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This paper describes a method for extracting rapidly varying, superimposed amplitude- and 
frequency-modulated signal components. The method is based upon the continuous wavelet trans- 
form (CWT) and uses a new wavelet which is a modification to the well-known Morlet wavelet to 
allow analysis at high resolution. In order to interpret the CWT of a signal correctly, an approx- 
imate analytic expression for the CWT of an oscillatory signal is examined via a stationary-phase 
approximation. This analysis is specialized for the new wavelet and the results are used to construct 
expressions for the amplitude and frequency modulations of the components in a signal from the 
transform of the signal. The method is tested on a representative, variable-frequency signal as an 
example before being applied to a function of interest in our subject area - a structural correlation 
function of a disordered material - which immediately reveals previously undetected features. 



PACS numbers: 61.43.Bn, 61.43.Fg, 02.30.Nw, 02.70.Hm 



I. INTRODUCTION 

Although the first example of a wavelet basis dates 
back to 1910 [Q, it was not until the early 1980's, with the 
work of Goupillaud, Morlet and Grossmann |Q| in seismic 
geophysics, that the wavelet transform (WT) became a 
popular tool for the analysis of signals with non-periodic 
characteristics, termed non-stationary signals |J, 

The WT allows a signal to be examined in both the 
time- and frequency-domains simultaneously. The WT 
as a time-frequency method has replaced the conven- 
tional Fourier transform (FT) in many practical appli- 
cations. The WT has been successfully applied in many 
areas of physics || including astrophysics, seismic geo- 
physics, turbulence and quantum mechanics, as well as 
many other fields including image processing, biological 
signal analysis, genomic DNA analysis, speech recogni- 
tion, computer graphics and multifractal analysis. 

The term WT is conventionally used to refer to a broad 
selection of transformation methods and algorithms. In 
all cases, the essence of a WT is to expand the input 
function in terms of oscillations which are localized in 
both time and frequency. 

Different applications of the WT have different require- 
ments. Image compression, for example, often uses the 
discrete WT (DWT) to transform data to a new, orthog- 
onal basis set where the data are hopefully presented in 
a more redundant form [Q. Other applications, particu- 
larly signal analysis ||, use the CWT, sacrificing orthog- 
onality for extra precision in the identification of features 
in a signal. 

The principal aim of this paper is to present the WT 
in a form well suited to the analysis of one-dimensional 
signals whose frequency components have rapidly vary- 
ing frequency and amplitude modulations. In order to 
achieve this aim, we introduce a new 'tunable', complex 
wavelet. This wavelet is based upon the well-known Mor- 
let wavelet Q but is better suited to high-resolution anal- 
ysis. The features of the WT using the proposed wavelet 



are understood through an asymptotic stationary-phase 
approximation to the integral expression of the WT spe- 
cialized to the new wavelet. We demonstrate the prop- 
erties of the WT using two example functions, a mathe- 
matical function and the other a realistic physical model 
function. 



We are interested in exploiting the complex WT in 
the analysis of structural correlation functions which de- 
scribe the atomic structure of disordered materials. As 
these correlation functions have different spatial regimes, 
they may be classed as non-stationary signals. Despite 
the overwhelming success of the WT in other fields, we 
are aware of only one other paper on this application of 
the WT. In that application, Ding et al. J(| studied an 
experimentally observed correlation function of vitreous 
silica using the Mexican Hat wavelet. We improve upon 
this single, prior application in three significant ways, 
namely the use of the complex WT, a 'tunable' wavelet 
and the method of interpreting the resulting transforms. 

We then analyze the reduced radial distribution func- 
tion (RRDF) of a structural model of a one-component 
glass with pronounced icosahedral local order The 
resulting WT clearly shows the existence of different fre- 
quency components in the RRDF and their exponential 
decay. These features were not clearly detectable by ear- 
lier methods (c.f. Ref. ||). 

In Sec. II we review the mathematical framework of 
the wavelet transform and discuss some mother wavelet 
functions before modify ing an existing wavelet for our 
purposes. In Sec. Ill A we consider the WT of a gen- 
eral oscillatory signal using the new wavelet. The re- 
sults are then used to interpret the wavelet tra nsform s of 
a variable-frequency examp le function in Sec. Ill B and 

Concluding remarks can be 



of the RRDF in Sec. |IIC 
found in Sec. [i^. 
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II. FORMULATION 

The underlying WT used in this paper can be com- 
pletely described as a one-dimensional complex, contin- 
uous WT using wavelets of constant shape || . We begin 
by examining the formulation of this WT in terms of an 
integral transform, before examining the choice of mother 
wavelet function. 

For simplicity we use time-frequency terminology, con- 
sidering the signal to be an input function of time, f(t). 
The CWT is an integral transformation which expands 
an input function /(f) in terms of a complete set of basis 
functions £(f; a, b). These basis functions are all the same 
shape as they are defined in terms of dilation by a and 
translation by b of a mother wavelet function ip(t): 

e(«;M) = iarv(^p). a) 

with a,b and a^O. 

The CWT F(a, b) is defined as the inner product: 

F(a, b) = (C\f) ee la]- 1 J°° r (^—^j /(f) dt. (2) 

The original formulation of the CWT used a pref- 
actor \a\~2 to give a normalization to unity, = 1). 
We choose an alternative prefactor | a. | 1 (giving = 
| ci | 1 ) after Delprat et al. [9, As we shall see, this for- 
mulation of the CWT allows for simple frequency identi- 
fication by examining the maxima in the modulus of the 
CWT with respect to the scale a. 

In order to understand the CWT, it is useful to relate it 
to the FT. The FT has a non-localized, plane-wave basis 
set and, therefore, has a single transform parameter - the 
frequency v. In contrast, the basis set of the CWT con- 
tains localized oscillations characterised by two transform 
parameters - the scale (or dilation) a and the translation 
(or position) b. It is this critical difference which makes 
the CWT preferable for the analysis of non-stationary 
signals. 

We are free to choose a functional form for tp(t), sub- 
ject to some constraints. Some of these constraints are 
forced upon us whereas others arise from the practical 
usefulness of the resulting transform. 

In order to recover a function from its wavelet trans- 
form via the resolution of the identity ip(t) must sat- 
isfy an admissibility condition. Although we do not make 
direct use of the resolution of the identity in this paper, 
we require that our choice of if)(t) satisfies this condition 
to ensure that all information about the signal is retained 
by the transform. The admissibility condition is essen- 
tially that the FT, ^(w) = (e^*]^), satisfies the relation 

■0(0) = 0, equivalent to requiring that the mother wavelet 
and, hence, the basis wavelets, have a mean of zero. 

Beyond simply satisfying the admissibility condition, 
it is practically useful to create mother wavelet functions 
which mimic features of interest in the signal. In the case 



of time-frequency analysis, mother wavelet functions are 
chosen which represent localized sinusoidal oscillations. 
The resulting wavelet transforms can then be used to 
extract instantaneous measures of frequency and ampli- 
tude. The uncertainty principle dictates that the prod- 
uct Af Acu of the time and frequency uncertainties of such 
wavelets has a lower bound. It is no surprise, therefore, 
that this class of mother wavelet functions are typically 
based upon Gaussians. However, it is still possible to 
trade temporal precision for frequency precision by al- 
tering the number of oscillations in the envelope of the 
mother wavelet. 

The simplest such wavelet is the "Mexican Hat" wavelet 
which mimics a single oscillation and is commonly used 
in signal analysis. The functional form of this wavelet 
is the second derivative of a Gaussian. This wavelet of- 
fers good localisation in the time domain whilst retaining 
admissibility. However, this wavelet has two major draw- 
backs for general signal analysis: (i) useful information 
can only be extracted from the WT at discrete intervals 
where the wavelets are in phase with the signal, and (ii) 
the time-frequency resolution is fixed. 

The former drawback has been overcome by the in- 
vention of complex wavelets which mimic localized plane 
waves. The WT can be computed separately for the 
real and imaginary parts, yielding a complex scalar field, 
F(a, b), where the modulus and argument of F represent 
the amplitude and phase of the signal, respectively. 

The latter drawback has been overcome by the inven- 
tion of tunable wavelets which include an additional pa- 
rameter to the mother wavelet function controlling the 
number of oscillations in the envelope. 

Goupillaud, Morlet and Grossman overcame these 
problems simultaneously with the invention of a mod- 
ulated Gaussian wavelet, now known as the "Morlet" 
wavelet This wavelet has a parameter, a, which con- 
trols the number of oscillations in the envelope, allowing 
time and frequency uncertainties to be traded. Thus the 
Morlet wavelet can be expressed as: 

ihi(t;tr) = n-ic M ((T)e-i t2 (e icTt - n{a)) , (3) 

where c M {<r) = (1 - 2e~i a2 K,(a) + k^o-))" 1 / 2 and the 
parameter n(a) allows the admissibility condition to be 
satisfied. 

The FT of this wavelet is: 

^■a) = Tr-ic M ((T)(e-i^ ^ - K (a) e~^) (A) 

From Eq. (0) it is clear that the admissibility condition 
■0m(O;ct) = implies that k(c) = e~2 a . 

Many previous applications of the Morlet wavelet have 
been concerned with signals containing slowly varying 
frequency and amplitude components for which large val- 
ues of a (> 5) are applicable and n(a) (< 10~ 6 ) is negli- 
gible §. 

However, we are interested in applying this type of 
analysis to signals which contain rapidly varying frequen- 
cies and amplitudes. In this case, the ability to use small 
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Figure 1: Envelope |?/<m| 2 of the Morlet wavelet ipM{t;a) 
(Eq. [i)) showing the unwanted transition from mono-modal 
to bimodal behaviour at small a (< 1.79785). 
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Figure 2: Morlet wavelet ipM(t;a) (Eq. []) for a = 1: real 
part, solid line; imaginary part, long dashed line; and envelope 
±|t/>|, short dashed lines. 



values of a becomes important as we wish to maximize 
the temporal resolution by minimizing a whilst still be- 
ing able to separate the various frequency components in 
the signal and, consequently, n(a) is no longer negligible. 

Although the Morlet wavelet is admissible at small 
ex, the temporal localization is unsatisfactory (see 
Fig. 0); namely, \iPm\ 2 undergoes a transition from mono- 
modality to bimodality (a single ridge at large a splits 
into two symmetric ridges for small a). The wavelet 
transform of a signal performed using a wavelet which 
has a bimodal envelope results in the signal being local- 
ized about two different positions (see Fig. ||). This pro- 
duces unwanted artefacts in the resulting instantaneous 
frequency and amplitude measurements (shown later in 
Figs. |l| and 0). 

Therefore we remedy this drawback by modifying the 
Morlet wavelet to produce a new wavelet, 4>{t;a), such 
that IV'I 2 has a single, global maximum for all a. For the 
new wavelet we choose to replace the single, normaliza- 
tion constant cm(c) in the Morlet wavelet with two new 
parameters p(a) and q(a) determined by two conditions: 
(i) total normalization of the wavelet to unity, and (ii) 
equal contributions to the normalization from the real 
and imaginary parts. The new wavelet has the following 




Figure 3: Envelope \ip\ 2 of the new wavelet (Eq. 

functional form: 

ip(t; a) = 7rX e X* 2 [p(a) (cos(erf) — /c(cr)) 

+ i q(a) sin(cri)] , (5) 

where p(a) and q(a) are given by: 

p(a) = (l-e V ) 4 , (6a) 

q{a) = (l + 3e- CT2 -4e-3 CT2 ) " . (6b) 

The Fourier transform of this wavelet is: 

fee) = i e -^+"> 2 (e™ - 1) 

xr-l)p(<7) + (e'Xl)9(a)]. (7) 

It is worthwhile noting that the real part of this new 
wavelet recovers the functional form of the Mexican Hat 
wavelet in the limit a — > 0: 

R*#(i;0)] = yf^e-5 t2 (t 2 - 1). 

Thus the new wavelet allows a complete transition from 
very high temporal localization, a — + (the Mexican 
Hat wavelet), to maximum frequency localization, a — > 
00 (plane wave). Even in the limit of minimal a, 2 
remains mono- modal (see Figs. || and ||). Thus we have 
improved upon the temporal localisation of the Morlet 
wavelet. 

We have also checked that, using the new wavelet, the 
original signal can be recovered by the resolution of the 
identity operator. 

III. ANALYSIS 
A. Instantaneous Frequency and Amplitude 

In this section we demonstrate how the new wavelet 
may be used to extract instantaneous frequencies and 
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Figure 4: New wavelet ip(t;a) (Eq. for a — 1: real part, 
solid line; imaginary part, long dashed line; and envelope 
ilV'L short dashed lines. 



amplitudes from a signal via the CWT. The following 
analysis is based upon the stationary-phase approach of 
Delprat et al. H but is specialized to the new wavelet. 

The wavelet transform F(a, b) at a given scale a and 
translation b is given by the integral (Eq. ||) of a rapidly 
oscillating integrand. This integral may be rewritten in 
the form: 



F(a,b) 



where: 



A(t; a, b) 
<S>{t;a,b) 



a i$(t:a,b)+\nA(t;a,b)^ ^ 



A f (t)At(t;a,b) (9a) 
(j>f{t)-^(t-a,b) (9b) 



with f(t) = Re[A f {t)e^f^] and £(t;a,b) = 
A^t;a 1 b)e l4 '^ a - b \ 

In order to take the integral in the stationary-phase 
approximation, we first approximate A% by a Gaussian. 
From Eq. H we have Aj(t; a, b) = a~ 2 A^((t-b)/a), where 
Ay is taken to be a normalized Gaussian whose variance 
(T^(ct) is equal to the variance of \if)\ 2 , giving: 



A 2 Jt;a,b)~\a\- 



27TC7, 



■ exp 



(10) 



where the variance a^(a) can be found analytically: 
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The approximate envelope, A 2 , tends to the true en- 
velope, |?/>| 2 (see Fig. ||). 

We assume (without loss of generality ||) that there 
is a single point of stationary phase for the integrand in 
Eq. (||) at t — t s (a,b). Under the conventional asymp- 
totic approximation: 



d<P f 



Of 



> 



1 dA 



A dt 



(12) 



we expand $(t) around the stationary point t s assuming 
<&"(t s ) ^ and substitute the approximate expression 
for A$ from Eq. (^) into the integral, which can then 
be taken. This gives an approximate expression for the 
squared modulus of the CWT using the new wavelet: 



\F(a,b)f 



-a^ 2 (^)(l + 4a 4 4ci>"(t s ) 2 ) 



x exp 



a 2 a 2 ^'(t s ) 2 (t s -b) 2 
l + 4a 4 ( j4$"(t s )2 



(13) 



Further, assuming the frequency of the mother wavelet 
to be constant {4>'i(t) — 0) and the frequency varia- 
tion of the signal to be slow in the region of interest 
(i.e. \$>"{t s )\a 2 a 2 < 1) then: 



\F{a,b)t 



a^A 2 (t s )e 



-aM*"(t s ) 2 (t s -b) 2 



(14) 



For a monochromatic signal (i.e. a signal which con- 
tains only a single frequency at any given position) , there 
is a scale a r (b) at any given b which corresponds to a ba- 
sis wavelet centred at b whose frequency <^(&; a r (b),b) is 
equal to the local frequency of f(t). The scale a r (b) of 
this wavelet identifies the instantaneous frequency of the 
signal and may be found as the solution of the equation 
$>'(b;a r (b),b) = 0. From the definition of the points of 
stationary phase ($'(t s (a, b); a, b) = 0), this corresponds 
to t s (a r (b),b) — b, an alternative equation which can 
be used to find a r (b). With the choice of normalization 
used in Eq. (|l|) , it is clear that these points maximize the 
expression for the squared modulus of the CWT with 
respect to a as obtained by the stationary-phase approx- 
imation, Eq. 

As the CWT is a linear operation, superimposed fre- 
quency components are manifested as different scales 
df\b) which locally maximize \F\ (assuming sufficiently 
large a to resolve the peaks). The curves formed by the 

points at? \b),b are known as the "ridges" of the trans- 
form The trajectory of each ridge can be used to 
extract the amplitude and frequency modulations of the 
corresponding signal components. 

An approximate expression for the instantaneous am- 
plitude, A f{t), of a signal component can be obtained 
by rewriting the stationary-phase approximation to the 
squared modulus of the WT (Eq. ^4() on the ridge, 



n K~j2 a^Af (b) 2 , in terms of A f (t): 



Af{t) 



1 2 

2™* 



F(aW(t),t) 



(15) 



There are two different well-known approximations to 
the instantaneous frequency <^(t)/27r. As each has rela- 
tive merits, we consider both. 

The simplest approximation to the instantaneous fre- 
quency is the rate of change of the phase of the CWT 
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Figure 5: Gaussian approximations (solid lines, Eq. ^) to the true envelope |?/>(t;<j)| 2 of the mother wavelet function (dashed 
lines, Eq. |) for: a) lim a -> 0; b) a = 2. 



with respect to b, evaluated at a r (b),b: 



« M 5 

Z7T 



^Arg^a® (&),&)] 



(16) 



The derivation for this expression using the new wavelet 
is identical to that of the Morlet wavelet given by Delprat 
et al. §]. 

The other approximation to the instantaneous fre- 
quency uses the equality of the frequency of the signal 
and of the wavelet on a ridge to create an expression for 
the frequency of the signal as a function of the scale a r (b) 
on the ridge and the frequency of the mother wavelet, 



(i) 



(t) 



OJtP 

"2V 



(17) 



Conventionally, ui^ is taken to be the underlying modu- 
lating frequency, a, of the mother wavelet function. How- 
ever, this is a poor approximation at small a. Therefore, 
the obvious definition of w$ is the modal average (posi- 
tion of the highest peak) in the Fourier power spectrum 

2 . Unfortunately, this expression for c*ty cannot be 
found analytically. However, even at small a, the spec- 
trum | ip\ 2 is nearly symmetric about the main peak (see 
Fig. ||). Therefore, the mean average is always a good ap- 
proximation to the modal average (see Fig. |?|) and, unlike 
the mode, the mean can be found analytically: 



u^a) = V¥ap(a)q(a) (l-e-3- 2 ) . (i 8 ) 

Using this expression for in conjunction with the 
relationship between scale and frequency in Eq. (17), a 
CWT may be plotted as a function \F(u)$ /2'kv, t) \ of time 
and frequency. 

Delprat et al. |J proposed that the phase-based instan- 
taneous frequency, Eq. (fl6|), is more accurate than the 
modulus-based measurement, Eq. (|r7|), and suggested 
an iterative algorithm for extracting signal components. 
Carmona et al. have since shown that the modulus-based 
measurement is extremely resiliant to noise Jin] and have 




Figure 6: Fourier power spectrum \t[>(io, a) | 2 of the mother 
wavelet function (c.f. Eq. |7|) in the limit a — > 0. 




Figure 7: Approximations to the frequency uty of the mother 
wavelet function (Eq. ^): Analytic mean of |?/>| 2 (solid line, 
Eq. ^|), numerical mode of jf/>j 2 (long dashed line) and the 
asymptotic approximation oty = a valid in the limit a — > oo 
(short dashed line). 



suggested numerous methods for extracting signal com- 
ponents using this approach lljj] . 

Thus the instantaneous frequencies and amplitudes of 
components in a signal may be found from the CWT 
at the points where \F(a, b)\ is locally maximized with 
respect to a. These maxima can be identified numerically 
from a set of samples of F(a,b) generated by discrete 
approximation to the integral expression for the CWT, 
Eq. (j|). Once found, the maxima may be interpreted 
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Figure 8: Variable-frequency function f(t) (Eq. |lli| ). 




Figure 9: Modulus of the CWT as a function of frequency 
F{u3^,/2irv,t) of the function f(t) (Eq. |l|) using the new 
mother wavelet function ip(t; a) (Eq. ^) with a — 2. 

using the approximate analytic results given above. 

B. Example Function 

The method described in the previous section is most 
easily clarified by the following examples. First, we 
choose to apply the method to the simple, variable- 
frequency function (see Fig. ||): 

f(t) = sin(i 2 ) . (19) 

The FT f(uj) = \ (cos(±w 2 ) - sm{\u?)) conveys little 
useful information about the original function. 

However, the modulus of the WT does convey useful 
information, particularly when plotted as a function of 
frequency instead of scale (see Fig. ||) as this highlights 
the linearly changing local frequency of f(t) (given by 
Vf = \d/dt 0/|/27r) as a function oit. 

The CWT of f(t) contains a single, 'V shaped ridge 
at a r (b). This ridge reflects both the frequency modula- 
tion of f(t) (see Fig. [l(]) and the amplitude modulation 
(see Fig. |ll| ). In all cases, the results show fluctuations 
linked with the phase <pf of the signal. However, com- 
pared to the Morlet wavelet, the new wavelet produces 
much smaller fluctuations in all results. 



C. Reduced Radial Distribution Function 

We now apply the method described in Sec. 3.1 to a 
function of practical interest. We choose to study the 
RRDF d(r) of a model glass structure. 

The RRDF analyzed in this paper is taken from a 
structural model of the icosahedral (IC) glass ||7| c reated 
in a classical molecular-dynamics simulation [|12|. We 
calculate the transform as detailed in Sec. 2 and perform 
the analysis as discussed in Sec. 3 in order to study the 
components of d(r). The function d(r) is considered to 
be zero outside the range < r < L/2, where L/2 ~ 25 
is half the side of the cubic simulation super-cell which 
contains 108,000 atoms. 

The RRDF, d(r), is defined in terms of the atomic 
density p(r) as: 

d(r) = 47rr (p(r) - p ) , (20) 

where pnis the average atomic density [jl3| . This is shown 
in Fig. yja for the IC glass. Reduced Lennard- Jones 
units (r.u.) are used for length with the mean nearest- 
neighbour separation being 1.15 ± 0.05r.u.. The damped 
extended-range density fluctuations are clearly visible, 
extending beyond lOr.u. (see the inset in Fig. |l2|a). 

From the Fourier power spectrum of d{r) (shown in 
Fig. |l2|b), it is clear that d(r) contains many components 
with different frequencies. The highest peak in |<i| 2 occurs 
at the frequency v d — 1.08 ± 0.01. This peak has non- 
zero width implying that the real-space fluctuation in 
d(r) corresponding to this peak has a spatially varying 
amplitude but we cannot deduce a functional form from 
this alone. 

Plotting the modulus \F d (a,b)\ of the CWT using dif- 
ferent envelope widths, shown as a function of r and 
v (= 27rfc) in Fig. [l3|, allows d{r) to be examined in 
the time-frequency plane. Using small a results in high 
spatial resolution but poor frequency resolution and the 
ridges are smeared together (see Fig. fl3|a). Larger val- 
ues of cr separate the ridges at the cost of decreasing 
the spatial resolution (see Fig. |l3|b). Unlike the exam- 
ple function from the previous section, d(r) contains sev- 
eral components with different frequencies which, partic- 
ularly when using large a, manifest themselves as sepa- 
rate ridges in the WT. In this paper we consider only the 
prominent ridge at v ~ Vd but the same analysis can be 
applied to the ridges seen at other frequencies. 

The ridge along v ~ Vd shows that the prominent 
frequency component identified from the Fourier power 
spectrum of d(r) is particularly strong around r = but 
decays away with greater r. This trajectory of the ridge 
can then be used to extract the instantaneous frequencies 
and amplitudes of this component in d{r) . 

The instantaneous frequency found using a = 3 (see 
Fig. |l4]) remains constant over a large range of r. As 
expected, the scale at which this ridge occurs in the CWT 
of d(r) corresponds to the position of the prominent peak 
in the Fourier power spectrum of d{r). 
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Figure 10: Instantaneous frequencies Vf{t) using the new wavelet (solid lines), Morlet wavelet (short dashed lines) and expected 
value v(t) = \t\/n for |t| S> (long dashed lines) of the example function f(t) (Eq. |l^) with a = 1. Extracted using: a) maximal 
|F|, Eq. ||); b) 8ATg[F]/db, Eq. ©. 




Figure 11: Instantaneous amplitude Af(t) (Eq. 09) of fit) 
(Eq. |l9| ) with cr = l. Instantaneous amplitude A~f{t), solid 
line; expected amplitude of 1.0 for \t\ 3> 0, dashed line. 



The amplitudes of components in an RRDF are ex- 
pected to tend to zero in the limit r — ► oo for a disordered 
material due to the absence of long-range order. The in- 
stantaneous amplitude Ad(f) of the dominant ridge ex- 
tracted from the CWT using the new wavelet, Eq. dH|), 
is shown plotted on a logarithmic scale in Fig. 



15. The 



amplitude is clearly seen to decay exponentially in the 
region 2 < r < 18. The reason for the exponential form 
of this decay (also observed by Ding et al. || for silica 
glass) is not yet known. 

The method used by Ding et al. j6| could not reproduce 
the frequency modulation of the damped density fluctu- 
ations in d{r) and their observed amplitude modulation 
contained only six points which were noted to decay ap- 
proximately exponentially. In comparison, our method 
reproduces true, instantaneous frequencies (analogous to 
frequencies obtained by Fourier analysis), showing the 
frequency modulation of individual density fluctuations 
in d(r), and the amplitude modulations of these variable- 
frequency components, as a continuum of points. This 
gives much more compelling evidence for the exponential 
decay first observed by Ding et al. ||. In addition, we 
can detect a significant deviation from the exponential 
decay of Ad(r) at large r (see Fig. |l5|). This may either 
be due to statistical noise from the finite nature of the 



simulation or due to the use of periodic boundary condi- 
tions in the model producing effective long-range order. 
The precise reason needs further investigation. 



IV. CONCLUSIONS 

We have identified the complex, continuous wavelet 
transform using wavelets of constant shape § as a 
method well suited to the time-frequency analysis of 
one-dimensional functions. For our target application, 
namely the analysis of functions with components which 
have rapidly varying frequency and amplitude modula- 
tions, we have illustrated an important shortcoming of 
the existing Morlet wavelet explained the origin of 
this shortcoming and proposed a new wavelet which over- 
comes the problem. In addition, we have specialized 
an existing method [^) for extracting instantaneous fre- 
quency and amplitude measurements from signals to the 
new wavelet. 

Two example functions have been analysed using the 
new wavelet and new method of analysis. The first, a 
simple variable-frequency function, illustrates the signif- 
icant improvement of the new wavelet over the Morlet 
wavelet and gives numerical evidence that our method 
of analysis is accurate. The second is a real- world ex- 
ample of a direct-space atomic correlation function of a 
glass which highlights the advantages of the method over 
the conventional Fourier transform and greatly improves 
upon the single, previous wavelet analysis of such a func- 
tion by Ding et al. ||. 

We have successfully used the WT to analyze the re- 
duced radial distribution function (RRDF) of a model 
glass and can immediately identify previously undetected 
features. The dominant component in the RRDF (the 
damped extended-range density fluctuations) has a pe- 
riod which rapidly settles to a constant value. Other 
components with different frequencies are present in the 
RRDF. These oscillations all have approximately expo- 
nentially decaying real-space amplitudes. 
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Figure 12: a) RRDF of the IC glass (inset showing a magnification of the damped density fluctuations); b) its Fourier power 
spectrum |d(fc)| 2 . 



v 1.5 




v 1.5 



0.5 




Figure 13: Modulus \F d (ui^, /2-kv, r)\ of the CWT of the RRDF d(r) plotted as a function of frequency v (see Fig. ^) using the 
new wavelet (Eq. 0) for: a) a — 2; b) a = 15. 
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Figure 14: Instantaneous frequency vAr) of the largest com- 
ponent of the RRDF d(r) (see Fig. |l2]b). Solid line is v d ex- 
tracted using Eq. (^) with a — 3, dashed line is the best-fit 
constant frequency u d = 1.08 over the range 5 < r < 18. 
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Figure 15: Instantaneous amplitude Ad(r) of the largest com- 
ponent of the RRDF d(r) (see Fig. |l^) plotted on a logarith- 
mic scale. Solid line is Adit) extracted using Eq. (JL5|) with 
(7 = 2, dashed line is the best-fit exponential decay Ad(r) — 
ae~ l3r over the range 5 < r < 18, where a = 31.8 ± 0.5 and 
/3 = 0.350 ±0.001. 



